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2) B^E«_hH:ftg|5^U**IJfflLTB*-BU-li:Mff1--i>o 

3) E<OTHAg|5*^flJfflLTD*-L-lD:5.0^F — L-lD 
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else 

c ^ ^k^r^» ^ rowupdate ^^■lf^LX'^mti> 
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BIT (b) nmmMU^im^x-^m-t^o ^<D^k. mi u) mmu^^mi 

[0 0 3 2] 

^t;, 0 7 (c) mmu^^. 0 7 (c) ^^m.^mntm7 (d) *:MEg:ft 
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^<oW:. El 7 (e) mm^»'S:m7 (f) fiS«IHftSB*^^v»jEifi--2.o Sl^ 
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-5.0 ftmtc«iiit«i5^^iii7 (f) ±m=MU»i:m-ox-^m-r^o 

[0 0 3 31 

vitzl^i^ LT, mmi-^ J: -? LT ^>^v^o S mx.i£, B(7)#hs . Jiff 

i±tf:^r6] (^^m^T'^tij) L7tffl:^icrsBLTafe?ijtcMai-'?)o 

[0 0 3 4] 

[0 0 3 5] 

Mm3-K#BB<, .^(c, US (a) ^mmm'&'Srms (a) ^rcofws 

izgftSI5^t08 (a) :&<oii^si5^^«o-cssfi-^o -e-com, 0 8 (a) :£co 

mm^»t:m8 (a) ^<7)*:i^=:ft^*ffioTM«r-t--Sc 

-^Iw. US (b) &<Dm$.m»i. 08 (b) ^(^**tEgfti:08 (b) :£ofi 

-?-<^m, 0 8 (c) ^mmu^^ms (c) *w5s«aHfta5^^^sev^M«fi-'& 

o Mi;> 0 8 (c) ^Ot4j^S|5^=Sr0 8 (c) :^(DmmmAtm8 (c) swia 
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-S>o J':J.TOMm3- K-C\ :£3gUc t#v^TV^-5,tfi±3^ > h^TT-^^o array 
a(k,n)i±. :Sfr^T?iJ5r**:)^-<§^f?iJ<0S^^*&tt1--|>Bd^ijT'*)^o ip (n) {± 

^BB?iJ-t?*^o LU^IS<7)-9-r;W-'^>coT;i'3'i;XA(±, #iS¥l 2-358 
2 3 2-f-5r#HS$tL7tv^o nb(±, L U^HI-^l^l-ifSI-^ 7'n y c^iS:^^ 

[0 0 3 81 

l'^<7)1^^y'ay 9iiZ^y^^X(OLU^mm^t>^t. ip ( i ) 75« i J: 
v^ti;, 1f?iJOi^fi5: i p (i) firetAtt#x^o -?-LT, u p d a t e -9- 
r)U-^>if^Z/iiiLX. ^fylJOSSf^^f9o LU^)»7i-f>up da t e^y)l^ 

LffOo -C-LT, ftm<7)JiS7*n y i^lcov^TJi, SU^L U^Mi: u p d a t e 
Sr js C: o -CiffliS * *?-T1- o 
[0 0 3 9] 

10 0 4 0] 

9o Ctiti, e-updatel Vic J: oTfi'd o -^-eo^, TIrI 
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1 e nJi, lo(0;^Wy K:;4?iaa-t-^rn y 5'(^)illi-eib-5o i s l±, *8!M-ti, 
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d f - u p d a t e(±. "/n y DSt>'F S:MSf-f Vfib^o rn 

«:ff/i^:7n 7 ^'(Dftl^KDfiB (nbase2) tLTltttL, 1 e n ^fr/wicit» L 
, :/D y ^<7)fti5!)i:ftf^<7)fiS i s 2 t i e 2 ^ff^wlcft^L, d f - u p d a 
t eT-D t F<7)Mif*^fv\ >'N''J T|^WSrfl5i-&o 
[0 0 4 2] 

if-9->^;W-^>' e-update2 iVf^ViUL. ^^'JTWimiM^o ±fitmW. 
IC, 1 e n. is. i e ^rft#L. y B i: H<?5Mir-'^-■f-> b h - u p d 
a t e ^DftFaiL, ^Of*, nbase2 1 5|c*i> T , 1 e n> i s 2, i e 2 ^r**!), 
Sl/b h-u p d a t eT'^S^^rV^ ^^'J T|BI^^]5l*o 
[0 0 4 3] 

sEt^. I- y K<7)#-f-;6n 60^, Te <7)Mff3j t LT, e-update3 

M(c. ^<^^, ^(7)A^^#x./c*SO«]^.^7cUM-tfc^6. 1 e n. i 

s. i e ^tfMLfz'lk. •*)-y")U-^> exchange <t -3 r^y(7)A^X#£ ^ 

[0 0 4 4] 

EI 1 1 (i, rn -7 B t rn y DCDMff-'f7;u--^>o^'ta3- K-ei(5.2)o 
rnyi5'BOMilft^i5V»-Cii, -fr;!—^- > b - u p d a t e l±. **0^f^ij 
fi5?iJa(k,n)tCTi5'-b^ L, ±121^^ t H Ci^i^^^o 1 e n , isK iel^ 

J:o-C^T?iJa<0-7'D-y ^'BcoMif?r1-^o CClt?. i s : i e i: v^t fB-^(i, fi^ 
[0 0 4 5] 

rn y j'DcoMiftcjsv^-Cfi. •9-y;P-^>d-u p d a t etcJ:oT 
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[0 0 4 61 

H 1 2 ii. y n y ^ C (h A<7)MSf-9-:7^>'i/--f->'(7)^mn-- KT^^o 

u/ ^ C(7)Mif-9-7>-^>' c - u p d a t eX'it. y^a y ^BtF(Dm% 
UJ:oTrn'y i^C^rSfrt-^o a(i:iof. is2:ie2) ;6^^rn y ^ 
C?:SL. a(l:iof. iof + lliof+blk) Ti^^rn y a 
(iof + i:iof+blk. is2:ie2) -A^y^u y Ir^-to 
[0 0 4 71 

^ AOlSff*9-r;i/-^> a - p u d a t e "Cii. rn y i5^B t D-crn 
y ^ Ai^mi-^o a(i:iof> is2:ie2) T^^^y^^n ^7 ^ A. a (1 : 
iof. iof + i:iof+blk) ^^^u u^i^B. a(iofH-i:iof 
+ blk. is2:ie2) -^^y^x^ y ^ DT$)^o 
[0 0 4 8] 

y ^ Gcomm-^'X^^-'f'y g - u p d a t e tci5V>T(i. mi^O-^fy;!/- 
^>t\s\m^^. ruy^m^'^mi^itm. I^TffiM^t-'Sr^i- l en. i s2. 
ie2. i o f ^t?:ft»L. rny^G^rny ^D^H-eMifi-'2)o a (i 
of + i:n. is2:ie2) ^^ya y ^ GT** a ( i o f + 1 : n. i 
of + i:iof4-blk) ^^yu y^H^ a ( i o f + 1 I i o f + b 1 k. 
i s 2 : i e 2) n 7 ^ D ^ o 
[0 0 4 9] 

rn-y ^EOS*U<7)M»f-^fr;U--?^>e-u p d a t e IX^it. E(Dnn^^ 

j^Pfo^f^^ ( i . i + 1 : b 1 k) -CHSfi-'Bo 
[0 0 5 0] 

rn y ^ Ecr>2mE(D'^^'^y)\^-^y e-update2 X^ii^ rn y ^ E 
ft^i3^J^J.«i<7)^f^^ h;i.s (i. 1 : i-l) tMn&^mk(Dm^^ V)i^s ( 

i + 1 : b 1 k, i) T*MffL. 'r^ y^E<DT^m=fm<7:>mmm^Mn&^<D 
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[0 0 5 1 ] 

Mia::!- KT'*;E>o 

'7 ^ EOSf^<7)Mff-9--:/>'^-^> e - u p d a t e 3T*li. rn y ^E<7) 
JiH:fttf?«J*^tft^*i'mo^iJ--<^ h;^ s (1 : i-l> i) i:tf^^h>'^s ( 

i. 1 : i-1) -eM»fL> rny i:^EcoWftsmJ^:Jimr<^^^tiM*fmr<^«-^^ 

[0 0 5 2 ] 

rn '7 ^D:S.C/F(?DMif-9-^*>'^-^> d f - u p d a t e tlj3v>T(i. yn y 
^<^iliS 1 e n:j!in 0 J: »9/h?v^^^, -/n D*;^ vM^-T^u ^ F (•9--7';i/- 
^><7)?IS5:coi s. i e oT^^^^tt-?)) Sr^n ^7 ^ e (JOg^^jl s (1 : i 
- K i ) iig^cotf^^ h;i/a ( i . is! i e) TMIffl-^o ^Cl-C. :7^n 
•7 i5^D*>^vMiFOS*fit}&^^a (1 : i s : i e) h ^^oTV^JgjOli, 

B;6^^:^7-t 7 h $ttr|%<^jA^ tLTV^^/::^, - i - 1 <7)S^fiHO 

v^Tft»-r^Ci::6^^ro^7^D*-i,VMi:F(Cov^TW»-r-i)^i:ti^^o 1 en 
751^^2 0J2JL±. 3 2mTO*>g-tlfi. lenK len25r^ftL. df-upd 

a t e ^Wf e^tco^^^/m mm<0^mi Ltz^k. Mt^. df-update^ 

[0 0 5 3] 

IwmtcjoV^TIi. bh-update(i. le ni)^\ 0 i: ^A^^^^t^^ l^m 
<7)?S®(C J:oTMif 1 e n:6^^2 0m±. 3 2 J^J.TOB#tCfi. lenl. le 
n2i^mL. ^(r>m<Of^^\Z{t^ lenl. 1 e n 2 ^gUt^^a L. bh-u 
p d a t e $:DfU^^mL> IWIBlOj^tZ ct -^?S»:^^TV\ jgtcb h-update* 

[0 0 5 4] 
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-e. SUN<7)Sr^tt@^>f VSUN Performance libraryO^fgtlJt 
[0 0 5 5] 

[HI 6-112 9(i. ^iftn- KO^?!LS^-7n-^^- h -e«L7t: ^) coT-*>l)o 
> t LTsharedge^^iJA (k. n) iXti-t^o >^ x ^ S 1 0 ^^(i. 

m^L. nb= (n+iblk-l)/iblk^i3:SL. i=l^l5:S-r^o ^^x-y^^S 1 1 -^li. 
i;6%b-UC^Lv^:*>575-^flSffi-^o ^^x y^S 1 1 <7)*iJ»f75^^Y E S O^^tcii 
. ;^x^y7^S 1 7tciitfo y^x-y^S 1 lO«»f;^i^'N0O:*'^(C(i. y^.x^yT'S 
1 2tlioV^T. nbase^(i-l)Xiblk=S:f5:^L. Xx y^y'Sl 3 "e. A (nbase+1: 
n> nbase+1 :nbase+n) co^^^y'u y ^J^ipliblktrTn y ^ L L. If >^ n y 

i:^SC/^TjE:^^f?iJg|5*«:Mifi-^o 'T^T-. ip (nbase+1 :nbase+iblk) MfO 

T;l-n^^VXA(i:^tc*mMA>i^^aiML7t#ltaiiStiK«$nTv^^o ;^xy 
14tIi5V^T(i. Hf'7*;U--^>exchgrowS:IIfC/(±5L. A (nbase+1 :nbase+ib Ik. 
1 :nbase) SB^J^O^f^^ h >'^<OAtt#x. ?:^f d o i^^?*?-^. j 5rnbase+l:*^^>nb 
ase+iblk^-e^fL^iirT. ip (j) > } ^mfzirt^A (j,l:nbase) tA (i 
p(j),l:nbase) iXUm^^o :^y^ yy^S 1 5^:^^. •9-7*;l'--^>update?:af ^ 
-e. ^liicorn y ^5:Hff1--i>o ^^x y -/S 1 6T*(i:. i=i+li: LT. :^r' y'fS 
1 ItCM^o :^T- y:7'Sl 7'e(i. nbase=(nb-l) X iblks^fS* L. Xxy^S 
IStlioV^T. A (nbase+1 :n, nbase+1 :n) (D^^dzmLX^ LU^)»=^ff^o i 
p (nbase+l:n-l) titf OAtL#x.(7)'|ff#:6^^AoTV>^o >^-r 'y:7*S 1 9 IC±5V>T 
. •iJ-r;W--^>exchgrow§:D?rFaSL> A (nbase+1 :n-l, 1 :nbase) SI545^0^f 
^;l/^0AtL#x.5r^f x y :7''S 2 0 tiiSV^T. ■^r;^--^>updateS:Pf ^'e 
Mifi-'&o >^xy:7*S2 lT^(i. -^f V^;l/-•^>exchgcol^Pf^'e^|J-^^ h;i/^A 
tL#x.^o 't't*?-^. j ?:n7&^^) lo^oi^?> LT'^^fkS^-^T. ip(j) > j ?:«7t-r 
A (l:n, j) , A (l:n, ip(j)) ^ Aj|x#;i>^o x y^S 2 2 Tii. ite^ij 
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[0 0 5 6 ] 

mi 7Rl/m 1 Sit. -y-yVlz-^^update^JO^^n-^^- htr^-^o 
;^T- y 3 0 t::j3V^Tti. ■9"7^;^"-^>b-update. d-update. c-update. a- 
update. g-update5:II?t/i±ii-t nbase. iblk. nmARZ/XUy Koffi^-C 
*'5>numthrd. I- y KcO#-§-nothrd^S:ttzgt"o -9-7*;^-^ >b-update^O? 

i5^^rjgifi-'i>o >^'r y7'S3 lT'(i. -fl-r;!^-^ >d-update5:0f 
T'n u/iJ^SrjSffi-^o v^^-r rs 3 2 (Ci^v^T. -^-r;!^-^ Vc-update?:O^^T* 
y ^^-^m-t^o ■7"S 3 3 -Cii. r;^-'5^>a>update^IIf^t:^r^ 

y^i:'^mir:ho :^'t y y'S 3 ■9-r;i/-'^>'g-update^lif y ^ 

[0 0 5 7] 

6(C:fev^Tfi. nothrdOfi|;6n 0>?.l^^y K;i-5*-=^^JBri~^o r- y 3 6 O 
*lJBff^^^N0Ot#^tc{i:, X-r •yT'S 3 StiittTo ;^t' ^y T'S 3 6 <7)fl|ff^^^Y E S 
O^t-^tlti, y •7°S 3 7 tcjoV^r. '9'r;^--^>e-updateS:Bf >y 
e ?:MffL. ^7-7^5 3 8 tciltro :^r' y^S 3 8tr{i. #>^1/ 'y KP^T^/smj 
Tmmi:M^o ^^r-y-f S3 9tCi5V^T(i:. df-update"e#>?. y K"e^ffl-r^ 
(is) . I^-.'S (ie) iik^X. SttvK1"o len=(nbase+nuinthrd-l)/numthrd 
. is=(nothrd-l) X len+1, ie=min(nbase, nothrdx len) HM^-t^o :7 n y :^ <7) 
l-^7Cg7fesRistart=nbase+l. .y ^ipilen=iblki:i-^o x ^ */ S 4 0 
v^Tfi. "^f r;^-^>df-update^B^^-cr^ y ^ f i^m-t^o 
[0 0 5 8] 

>^-r'yyS4 mi3V>Tii. I^i:<*&,^r. *l-.^.^it^ LSltiSi-o nbase2=nba 
se+iblk, len=(n-nbase2+numthrd-l)/numthrd^ is2=nbase2+(nothrd-l) X len+1^ 
ie2=min(n,nbase2+nothrdX len) HMWir^o v^T" y 4 2 tCioi/^T Ji. "/n 
1^7C§5feM i start =nbase+K ^7 n ^ ^iliiien=iblki:i-^o -y-y^l^--^^ 
df-update?rOf^t:^ y f SrMSf1-'i>o ;^f- y T'S 4 3 T-fi. #;;^lx >yK 
WT-/s'UT|^MS:^^o ;^-r yy'^S 4 4nothrd<7)fi|:6no^l/ y K75^5;4-^*lJBf 
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4 4 cO*lJSlf;6^^YE SO^i-g-tcti, ;^ :7' S 4 5 tcjs v^T . >e- 
update2SrlI^^try^ y ^ e SrMif L. Xr* y y S 4 6 (Ciltfo Xf- y-7°S4 6 
tzi5V^T(i:. #X W -;/ KF^r*/N*'; TIWIW^^^o y T'S 4 7 jSV^TJi. bh 
-updateO#Xlx-;/ KtC$tt^-r*&.*:. I^.^.^rft^-f -5)0 len=(nbase+ 
numthrd-l)/numthed. is=(nothrd-l) X len+1^ ie=min(nbase, nothrdx len) ^tt 
n^^o y iJ^o i^7cg(7)56Kistart=nbase+K y ^ipSlen=iblkt 1"^ 
o ■9-:/;^-^>bh-update*lI?AT% y ^ b ^'^mir^ho :^y^y'fS 4 S{z 
±5V^-r(i:. bh-updateO#X 1/ 7 KtcStt^gl-^,'^.. l^.'S?:lt»i"'i)o 
[0 0 5 9] 

^"r y ys 4 9 t-V^T(i^ nbase2=nbase+iblk^ len=(n-nbase2+numthrd-l)/nu 
mthrd^ is2=nbase2+(nothrd-l) X len+l> ie2=min(n, nbase2+nothrdXnothrdX le 
n)?:fft^L. y :5^<7) l^7c§O^Siistart=nbase+l. y ^ipilen=iblkt 
1--i)o X-r^yS 5 Otri5V^T(i. bh-updateO^X :y KtcSttiSI-^&.'S. 
.^.^nn-t^o X-r-yySS ltci5V>Tii. ^Xlx^y Vmx-y^') rmm^S: t :ho 
XT-^yyS 5 2-e{i:. nothrd<7)fiI;^)UOXl^^7 K7^-5;&^^*lJIBf-r^o X-r^y^S 
5 2<7)flff7&^^NOOJ#>^tC(i. X-r-yT'S 5 4tim<^. WBf^^^Y E S <7)^&'g"tz{i 
. XT-'y-T'S 5 StCi^v^T. •^r;^-•^>'e-update3^IIf^t:^ ra -y^e^M 
ifft. X-r y yS 5 4 tlittfo Xxy ^'^S 5 4 UioV^Tti. ^X P ^7 KFslTv^U 

[0 0 6 0] 

HI 9(i. •9-7';i^--^>b-update<hd-updateC0jjniS4-^-t:7n-^'^'- h-T^ 
-5) o 

•9-r;u-'^>b-updatetfi±. Xx-y7'S 6 0H:fov^-c. #x 1/ -y VX^j^Uir^ 
l>k7c:E<Oithj^.. 5: ft*1" -2)0 len=(nbase+numthrd-l)/nuinthrd 
. isl=(nothrd-l) X len+1, iel=min(nbase, nothrdX len) ^ iof=nbase?:It^'t" 
Xx y 6 1 (liJV^Tti. A (isl:iel, iof+l:iof+iblk) =A (isl:iel, 
iof+l:iof+iblk)xTRU-U (A (iof+1: iof+iblk, iof+1: iof+iblk)) -^irifM-t^ 



tbSE#2 003-3067303 



#il 2003-074548 



19/ 



[0 0 6 1 ] 

^j-r;l.--^>d-update-e(i:. >^ x y :7'S 6 5 t^ioV^T. y KT53^ffli--& 

$p.^,. i^.^.itfW^-t^o 1"^t>-^. len=(nbase+numthrd-l)/numthrd. isl=(noth 
rd-l)Xlen+l> iel=min(nbase, nothrdX len) > iof=nbase$:lt#i"'2)o ^"ryf 
S6 6Tii. A (iof+l:iof+iblk, isl:iel)=TRL (A (iof+1: iof+iblk, iof+1: io 
f+iblk)) -IxA (iof+1: iof+iblk, isl:iel)^tt^i"'i>o TRL(iiE::^tT?iJ<30T= 

[0 0 6 2] 

[12 0 (i. •9-'7\»U-^>c-update<:0 7n-^A'- hT^^o 

i"^o nbase2=nbase+iblk^ len= (n-nbase2+numthrd-l)/numthrd^ is2 

=nbase2+(nothrd-l) X len+1^ ie2=min(n,nbase2+nothrdX len) ^ iof=nbase^:ft 
Mir^o y^f-y-fSl ixa. A (l:iof. isl:ie2) =A (l:iof, iof+1 :iof+i lb 
k) XA (iof+1 : iof+i Ibk, is2:ie2) ^ft^L. '^r)l^-^>i:m^f^o 
[0 0 6 3] 

mz l(i. •9-r;l/-^>'a-update^7)*7n-^-v- h-e*^o 

i'-So len=(nbase+numthrd-l)/numthrd> is2=(nothrd-l) X len+K ie 

2=min(nbase,nothrdXlen) . iof=nbase4:ft*i"^o y ^ S 7 6 tCjSV^T (i 
. A (l:iof, is2:ie2) =A (1: iof, is2: ie2) -A (1: iof, iof+1: iof+iblk) XA 
(iof+1 : iof+iblk, is2:ie2) ^nfWLX^ •^':r)U-^>^mf^o 
[0 0 6 4] 

02 2 fi. ■^r;^-^>g-update<^-7n-"^^- hX^^o 

W't^o ir^t>^^ len=(nbase+numthrd-l)/urathrd. is2=(nothrd-l) X len+K i 
e2=(nbase.nothrdxlen) > iof=nbase?:tt^i-'i>o ^^-r 7 ^ S 8 1 tCi5V^T(i. 
A (iof+l:n, is2:ie2) =A (iof+l:n. is2:ie2) -A (iof+1 :n, iof+1: iof+iblk 
) XA (iof+1: iof+iblk, is2:ie2) itfWL. ^y^^l^-'^y itk^f 
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[0 0 6 5] 

112 3 (i. •^y^)U-^>e-update(Dy\J-^-^- h*tr^^o 

8 7tci3V^T. S(l:i-1, i+l:iblk) =S (1: i-1, i+1: iblk) -S (l:i-l,i) XS (i, 
i+l:iblk) . i=i+lSrft#L. ^J-^;!/-^ > Srlfett-So 
[0 0 6 6] 

112 4 (i. •9-y;U-f->e2-update<5C> >^n-f-^- h'e^).i>o 

f-f, rn 7^Eco;feJi|5S<7)g^7i^^S (k. *) oft^oog^t J: d 

S5:-eS:tt'i>o X-r -7 ys 9 0 -e(±. i=ltf5^L. ^7 9 1 tii^v^T. i>i 
blk;i-S;i>-^flJif-r'i>o -7^5 9 1 OfiJifT^^^Y E S 

. tmp=1.0/S(i, i). S(i, i:i-l)=tmpXS(i, l:i-l). S(i+1 : iblk, 1 : i-l)=S(i+l: ib 
Ik, l: i-l)-S(i, l: i-1) xS(i+l: iblk, i). S(i+1 : iblk, i)=-A (i+1, iblk, i) xtmp 
. A (i,i) =tmp. i=i+l?r?S^LT. t" ^> T^S 9 1 HM'i)o 
[0 0 6 7] 

tt\ y ^E(D&±m<D^m*^s (k. *) <Dm^j}<Dmmt^j:^xdi>z^\ 

mT-^tt^o ;^T-^y :7-S 9 5T-(i. i=li:IS^L. :^Tyy''S 9 6X\ i>iblk;6- 
S5ii-4-flif-t'5)o ^^-ry T'S 9 6<7)fiJ»f:6^^YE ScOJt'&-(::(i. '^r)\^-^y^ik 
If^o >^T- 7 9 6<7)*lJgff50^^NO<7)^^ii(i. TsT^yy'S 9 7 S(l: 
i-l,l:i-l)=S(l:i-l,l:i-l)-S(l:i-l, i) xs(i, 1: i-1). S(l:i-1, i)=S(l: i-1) xS 
(i,i). i=i+l4'«#L. ::^f' yys 9 eiiZM^o 
[0 0 6 8] 

02 6(i. ^j-'7Vl/-"^>df-update<7)-7n--^^- h-e^b^o 

miE# 2003-3067303 



#SI 2003-074548 



21/ 



§<7)^llistart^yn ^7 ^tpslen-trStt^o y^S 1 0 O^V-lt. len<10;^^S 
:^^^mm-t^o y V'S 1 0 0(Dmm:A^l^O<om^lZ{t^ y >^^S l O 4 tl 
ittfo XT- 1 0 0<7)fl»f7i^^YE SCO^&^tCii. Xx -y T^S 1 0 1 tCiSV^T 
. i=lfcf$^L. X-r*yyS 1 0 2 tcj^v^T. i<=len:^^S;i^^WSff1'^o Xt--^ 

^^s 1 0 2 (Dmm^'i^ O(om^i^zii^ ^y^)v-^y^mi>ho x^^y ys 1 0 2 

0*lJ||f:e>^^YE S<7)*if^tZ{±, X-r -y V^S 1 0 3 tl^v^T. js=istart. je=istart 
-l+len-l> A (js: je, is, ie)=A (js: je, is, ie)-S (js:je, i) XA (istart+i-1 
,is:ie) ^JS»L. X x ^y ^ S 1 0 2 (CM'i>o 
to 0 6 9] 

Xt" y 7'S 1 0 4r*li. len>=32*^v>iilen<=20:^-5;i-=^*iJBf-r^o Xr-y/ 
S 1 0 4 OflJif:6^^N0O^#^tZ(i. Xx y^S 1 0 6T\ lenl=len/3. Ien2=len 
-lenl t LT. Xx y :7''S 1 0 7 [z^t^o XT' y 1 0 4 O^lJif^^^^Y E S co* 
>g^tc{i. X-r y T'S 1 0 5 tcfcv^r. lenl=len/2. Ien2=len-lenli: LT. X-r 
•yyS 1 0 7tCittro XT--y:7*S 1 0 7 X-\t^ •^•r;^--^ >bf -update 
tiB^r/tti-r (rn ^y ^co^sgoistartiiyn y ^ificoienl ^S:tti£1-) o x-r 
y'/SlOS'Cfi^ js=istartx je=istart+lenl-l^ js2=js-nbase^ je2=js2+len 
1-K js3=js2+lenl> je3=js3+len2-l> js4=istart+lenl^ je4=js4+len2-K A ( 
js:je, is:ie) =A (js: je, is: ie) -S (js2: je2, js3: je3) xA (js4, je4: len, is 
:ie) . istart2=istart+lenl$:a#LT. X t- y y s i o 9 tCiiOo XT -y^'S 

art2trn^y ^itiS<7)ien2^S:ttffi1-) > > ^*?-Tt-'i>o 

[0 0 7 01 

02 7(i. •9-y;^-'^>bh-update<7)'7n-^-=^'- ht?^>2>o 
CO-9--r;^-^>(iS'jfW-/n^*'T AT'*);^o i-f^ ^Xl^y KT-^JiSi-^® 

H^^i-^.if.:5.0^^$l'.'S^is. ie-r^Stt^o i/c. rn «y ^ E (7):feJiPS(7)S^;6^^S 

tartt rn y i5^ipS^len^$frt>2)o 
[0 0 7 1 ] 
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;^xy7'*S 1 1 Sr-it. leTKlO^^^^^irmWr-t^o ^^T-^y/S 1 1 5 <7)ffl»f:i^^ 
NOOtir^tcfi. y^x^ys 1 1 9tcBtfo V^S 1 1 5 Oflglf^^^^Y E S O 

^^-r yV'S 1 1 6 joV^r. i=l<!:iS:^L. Xx ^ys 1 1 7tcio\/^ 
i<=len:*-S:*^=S:*lJKff-r^o >^t- 7>^S 1 1 7 0*IJ»f?^^^N OO^^Hii. ^^7* 
;l'--^>S:att^o ;^T-'7 7*S 1 1 7mm:^^YES(Dm'^i^^i±. Xt^ y^S 1 
18tCi5V^T. j = istart+i-K j2=j-nbase> je=istart+len-l. je2=je-nbase. 
A (is:ie,j) =-A (is:ie,j) XS (j2,j2) . A (is:ie,j) =A (is:ie,j) -A 
(is:ie, j+l:je) XS (j2+l: je2, j2) ^rS^LT. Xx^T^Sl 1 7 H^tfo 
[0 0 7 2] 

X-r -y^S 1 1 9-e(i. len>=32^>^VMi:len<=2074-S;d-^*iJBfi-^o T^'r y 
S 1 1 QOfflBfTi^^NOO^^tKi:. X-r^y^S 1 2 1 *iJ»f:i^^Y E S co^ 

X-r yS 1 2 0 tliitTo X-r y7'S 1 2 0 tCi^v^Tti. lenl=len/2. 
len2=len-lenl^lt*L. X-r 7 y^S 1 2 2 tijltfo Xf" y 1 2 1 iliSV^T 
fi. lenl=len/3. Ien2=len-lenl SttS^ L> Xr-yT'S 1 2 2 tciitro T^'r vf 

s 1 2 2tci3v^r(i. ■9--r;w-'f->bh-update?:Wfe^uBfr;^^mL. i^mra y 

^O^Mistartt -7^0 -7 ^ipglenlSrS^tiSi-o Xt^^^S 1 2 3 ^Cii. js=istar 
je=istart+lenl-l^ js2=js-nbase^ je2=js2+lenl-K js3=istart+lenl^ je3= 
js2+len2-l^ js4=js3-nbase. je4=js4+len2-l. A (is: ie, js: je) =A (isiie, j 
s:je) -A (is: ie, js3: je3) xS (js4: jse, js2: je2) . istart2=istart+lenl 
mMir:ho XT- ^yS 1 2 4tcib^v^Tii. -9-r;l/-^>bh-update^S«Wtcllf 

[0 0 7 3] 

1112 8 (i. ■9-r;^--^>exchgrowO'7n-^^- hx^^^o 
Xx yy^'S 1 3 OtcioV^T. /^VrmmiM^o XT- 1 3 1 tliBV^T. 
#Xlx ;/ K-e^jai-^$&,6;. ^^-.i^^^ft^i-^o len=(nbase+numthrd-l 
)/numthrd^ is=(nothrd-l) X len+1^ ie=min(nbase, nothrdX len) . j=nbase+l5' 

tfUir^o X-r yy^S 1 3 2-r«±. j>min(n-l,nbase+iblks) •V^^^^^:^^im 
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(i, Ts-f yy'S 1 3 Si^zmtso ^"ry y'S I 3 3-X:lt. ip(j)>j T^^STi^^^AWt" 
-i>o XT^ y 1 3 3cO*iJKf:^^^NO<7)B#ii:. >^7^ yy^S I 3 5 tiittro ^^'r y'f 
S 1 3 3<7)fIJ»f7J>^^YE S<7)Ptii, X-r T'S 1 3 4 ICiJV^-C. A (j,is:ie) i: 
A (ip(j), is:ie) 4-AtL#x.. J^T^yT'S 1 3 5 tliltro y y S 1 3 5 
j=j+l^it^L. Xt^ y ys 1 3 2 \iZm^o :^'r y'fS 1 3 

[0 0 7 41 

Xt^ y 1 4 Otl±5V^T. /^mj T|5l*8^:?t^o X-r ^y 1 4 1 tliJV^T. 
#Xl/ y K-C^fi-r^^&.ir, *l'.«S^tt^i-^o i-'t*^"^. len=(n+numthrd-l)/nu 
mthrd> is=(nothrd-l) X len+K ie=min(nbase, nothrdX len) ^ j=n-l?:ttXL^ 
y y^^S 1 4 2 tiiitfo y^S 1 4 2-trii:. j<l^ffl»fi-'i)o Xf-y 

1 4 2 60«»f:6^^YESc7):^'g-tI{i, y^T-^/S 1 4 6 tcjltro ;^T' y'rS 1 4 2 
mmi)^^0(r>m^K{i^ ;^7^>y7'S l 4 3CiDV^T. jp(j)>j^fiJKfi-'2)o ^ 
T-^y-y'^S 1 4 3mm:^^^0(7)m^iiZl±^ XT^-y >^'S 1 4 5 tciltro ;=^f-y y^S 

I ^ 3<omm^^YES<Dm^^Z{t^ ;^-r yT^S l 4 4tci5V^T. A (is:ie,j) 
i:A (is:ie, jp(j)) S:AtL#x.T. ;^-ry7'S 1 4 5 (littfo X-r-yy^S 1 4 

5-e(i. ]=J-lt Lr. XT- y 1 4 2 tcM^o Xt" y :7'S 1 4 6 "^^(i. >'^U 

?y Its -SSS^ T ;^ V X A M L T (i: . i^T <7)iS:#» # fis $ tt v ^ ^ 
[0 0 7 5] 

G. H. Golub and C. F. Van Loan "Matrix Computations" The Johns Hopkins 
University Press, Third edition 1996 

2003-3067303 
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■y-ft. 

[0 0 7 6] 
[0 0 7 7] 

(ittiS) ±ie:^-S(±. Gauss-Jorda n O;;^-^^ -/o -y ^^l,Z3t 
[0 0 7 81 

itfr^iJ =S: ^ § fi^ 1*1 {c , ^^<^iE:^ff^ -ru y^-^ }t i- ^ ;^ x y 7' i: . 
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[0 0 7 9] 

#K 5 icf ett<7? y n V A o 
[0 0 8 01 

(#Se7) ±fe:^fe{i, Gauss-Jorda n<D:fj^^-fa y ^^iz^ 

(#fe8) ffirfe#7'n yi?s:^s!iL.-c3fe?ijst«:-r-5.$v><^:s-iij<7)itii{±, mn 

^ $ tt ^ c: i: 5r i: -t ^ #S£ 5 ie«cD 7° a V A o 
[0 0 8 1] 
[f&?^0?**l 

[mi] 
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[1112] 
[Si 3] 
[0 4] 
[05] 

[He] 

[117] 
[US] 
[1219] 
[HI 0] 

[0111 

2^i&?go*ifef^E<7)^ia='- K (-e-03) T-^-So 

[012] 
[013] 
[014] 
[015] 
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[016] 

[mi 7] 

[HI 8] 

Kco^S^r :7n->^'^ 

[mi 9] 

Mian — Kco^aSr*7n — -^-v 
[E12 0] 

[[112 1] 

[112 21 

[1112 3] 

m2 4] 

[13 2 5] 
Mian - K<7)^a5:"7n---f'^ 
[02 6] 

[02 7] 

Mian - h^<D^m^ytl-^^^ 

[02 8] 

Mian - K(7)^a^"7n-^-v 

[02 9] 
Mian - KoMa^:7n-^^ 



-hT-*L/::0 (^0 4) -e*>i>, 



- h-C^L/::0 (^(7) 1 0) Xi^^. 



- hX-mttzm (^O 1 2) T**/5), 



- htrSLfcm {^(D 1 3) -e*-?), 



-hT^L7t0 (-?-0 14) T*-i>. 
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1 0- 1 - 1 0-n 

1 1- 1-1 1 - n 

1 2 m^^^m 

1 3- 1 - 1 3-n 
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1/ 



mm 



11-1 11-2 



11-n 



J 2 



^3- 



1=1 



13-2 13-n 



T ^ 

10-1 10-2 



10-n 
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0 0 3 - 0 7 4 5 4 8 ^-i^ '. 2/ 

[SI 2] 



A 


6 


C 




K 




D 


F 








G 


H 


I 
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E 





a 














s;- b 



E2 



IMA] 



E 



E3 





a 












d 
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me] 



5/ 





/ 








B 






C 








D 






E 






F 








C 








H 






I 
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[US] 




ailiE# 2003-3067303 



[119] 
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-^-v: 8/ 



subroutine inversion ( a, k, n) 
shared array a(k,n),ip(n) 
create threads 
set r.o-hrd and numthrd 
c nothrd It::^ U K<0## 1 TH, numthrQ=ttTH U y K^OSg^gc) 

iblk«-:/n -y9<^ 
nb» (n+blk-l; /blk 
dc i«l,nb-l 
ntaase=(nb-1 ; *blk 

call LU (a, k, n,nbase, blk/ ip, nothrd, numthrd) 

! ^f^J^Mff-fa^Jll^^r^^lJtCfT^c ip(nbase+l :nbase-blk; 

! 2-3 5 8 2 3 2-^^ciH5di^i^-CV^5^o-z?^i^^^^J^. 

do i-nbase ' 1, nbase+blk 
if f ip (i) >i) ther. 

exchange {a (i, 1 :nbase) , a (ip (i) , 1 :nbase) ) 

endif 

enddo 

call update (a, Ic,n, nbase, blk, nothrd, numthrd) 
encdo 

noase=(nb-J ) *blk 
bl klast=n-nba s e 

call LU fa, k, n, nbase, blklast, nothrd, nunthrd^ 
call update (a, :<,n, nbase,blklast, nothrd, nujithrd) 
return 
end 
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[mo] 



subrourine update (a, n, nbase, r>lk, notlird, numrhrd) 
call b-update (a, :<,n, nbase, blk, nothrd,marathrd) 
call d-update{a, :<,n, nbase,blk,nothrd,nunthrd) 
call c-update (a, !<, n, nbase, blk, nothrc, nunthrd) 
call a-update (a, <, n, nbase, bik, nothrc, nunthrd) 
call g-update (a, n, nbase, blk, nothrd, nuir.thrd) 
BARRIER SYNC 
if (nothread=l> then 
c e (D-^gr 1 

call e-updatel (a (nbaserl^nbase+l ) , k, ri,blk) 
endlf 

BARRIER SYNC 

ler.«— { nbase+nunthrd- 1 ) /numthrd 
is*— (nothrd-1 ] *len+l 
nothrd^len 

call df-update (a (nbase+1, !},:<, n, a {nbase-l, nbase-1) , is, ie,blk) 

nbase2'— nbase -bl k 

len*— (n-nbase2-numthrd-l ) /numthrd 

is2*-nbase2- r.norhrd-l) *len+l 

ie2*— nbase2-nothrd''ler. 

call df~update (a (nbase+1, 1} , k, n, a (nbase+1, nbase- 1 ) , is, ie, blk) 
BARRIER SyN:j 
if (noihread=li then 
c e(^H^2 

call e-upGate2 (a (nbase+1, nbase* 1 i , k, n, blkj 
endif 

BARRIER SYNC 

len*— (nbase+ numthrd- 1 J /niimtnrd 
is'^ {noihrd-i:» *len+l 
ie*~iicthrd*ien 

call bh-update {a (l,nba3e+l) , k, n, a {nbase+1, nbase+1 ) , is, ie,blk) 
nbasc2*— nbase+blk 

len*^ (n-nbase2+nuir.zhrd-l ) /nuitthrd 
is2*-nbase2+ (nothrd-: ; *len-l 
ie2*— nbase2+nothrd*len 

call bh-update {a ;i,nba3e+l) , k, n, a Jnbase-rl, nbase+1 ^ , is, le,blk; 
BARRIER SYNC 
i f (nothread-1 ) then 
c € COMICS 

call e-update3 (a (nbase+1, nbase+1 ) , k, n, blkj 
endif 

BARRIER SYNC 

leii= Jn+namthrd-1 ) /nixathrd 
is= :ncthrd-l) *len+l 
ie=:nin (n, ncthrd*len) 
if (ip (i) >i) then 

exchange {a (is :ie, i) , a (is, ie, ip(i) i • 

endif 

enddo 

BARRIER SYNC 

eliminate threads 

return 
end 
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subroutine c-upda-e ( a, n, nbase, blk, ncthrd, numwhrd) 

shared array a[k,n) 

len*— (nbase+numthrd-1 } /r.mr.thrd 

isl^ <nothrd-l) *len 

iel<— nothrd*len 

icf=ribase 

a (isrie, iof +1 :iof +blk) 

a(is:ie, iof-1 : iof +blic} *TRU-U (iof +1: icf +blk, iof -1 : iof +blk} 

return 
end 

c -f^ D COSEff 

subroutine d-updete {a, k, n, nbase« blk, nothrd, r.umthrc) 

shared array a{k,n) 

ler.*— (nbase + numthrd-i; /numthrd 

is 1«— (nothrd- 1 ; *len 

iel<-nothrd*ienisl 

iof=rbase 

G (is :ie, iof+1 :iof H-blk) 

♦-TRL (iof +1 : icf+blk, iof-1: iof iblk) "^*a {is : ie, iof-i : iof — blk) 

I TRL i-^r^nnv^i 

return 
end 
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c c o)^Sfi 

subroutine c-update (a, k, n, nbase, blk, nothrd, numthrd) 
shared array a{k,n) 
nbase2*~nbase+blk 
iof =nbase 

len*— (n-nbase2-fnumthrd-l) /numthrd 
is2*-nbase2+ (nothrd-1) *len+l 
ie2<— nbase2+nothrd* len 
a{l:iof ,is2:ie2) 
•-a{l:iof , is2 :ie2) 

-ail: iof , iof +1 : iof +blk) *a (iof +1 : iof +blk, is2 : ie2 ) 

return 

end 

c :/d 5/^ A cOMiff 

subroutine a-update (a, Vl, n, nba£»e, blk, nothrd, nuinthrd) 

shared array a(k,n) 

len*— (nbase+numthrd-l ) /numthrd 

132-^ (nothrd-1) *len 

ie2*— nothrd* len 

iof =nbase 

a (1 : iof, is 2 : ie2 ) 

^a(l:iof, is2:ie2) 

~a ( 1 : iof, iof +1 : iof +blk) *a {iof +1 : iof +blk, is2 : ie2) 

return 

end 



ffiiiE# 2003-3067303 



mis] 
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C ^J^D y ^ G t7>3e^ 

subroutine g-update (a, k, n, nbase, blk, nothrd,nuinthrd) 

shared array a(k,n) 

ien*— (nbase + r.umthrd-1) /numthrd 

is2*-(nothrd-l) *len 

ie2*— nothrd* len 

iof^nbase+blk 

a (iof +1: n, is2 :ie2) 

*~a (iof+l:n, is2: ie2) 

-a (iof +1 : n, iof +1 : iof +blk) *a ( iof+1 : iof +blk, is2 : ie2) 

return 

end 

subroutine e-updatel (s, k,n,blk; 
shared array s(k,n) 

do 1=1, blk 

s (l:i-l,i+l :blk)*-s {1:1-1, H-l:blk)-s(l:i-l,i)*s(I,i+l: blk) 

enddo 

return 

end 

c 'T'a *y ^ E O 2[sl@COM^ 

subroutine e-update2 (s, k, n, blk) 
shared array 5(k,n) 

do i=l,blk 

tmp*-l .0/s (i,i) 

s (i, l:i-l)*-tmp*s (i, 1 :i-l) 

s(i+l:blk, l:i-l)— s(i+l: blk, l:i-l)-s (i, 1:1-1 )*s(i+l:blk, i) 

s(i+l:blk, i) « — s (i+1 :blk,i) *tmp 

s (i, i) *~tnip 

enddo 

return 

end 
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subri-urine e-update3 (s, k, r., blk) 
shared array s(:<,n) 

do i=l,blk 

s s l:i-l)-s(l:i-l, i) *s (i, 

enddo 

return 

end 

subroutine df~ update (a, k, n, s, is, ie, len) 
shared array a ( k, " ) , s { k, * > 
if (len<10) then 
do i~l,ien 

a (1 : is : ie) •»— a (1: is : ie) -s (1 : i-1, i) *a (i, is :ie) 

enddo 

else 

ifaen>=32 or ler.<=20 ) then 

ienl*~len./2 

len2*— len-lenl 

else 

lenl*— len/3 
ien2*— len-ienl 
endif 

call df-update (a, k,n, s, is, ie, lenl) 
a (i : lenl, is : ie) 

*— a (1 : lenl, is: ie) -s ( 1 : lenl : len 1+2 :len] *a (lenl+1 : len, is, ie) 

call df-update {a{lenl+l,l),k,n,s (lenl + 1, lenl+2 ) , is, ie, len2) 

endd f 

return 

end 
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c B 4oJ:t)^ H 

dubroutine bh-update (a, k.,n, s, is, ie, len) 
shared array a ( k, * ) , s ( ic, *) 
if ( len<10} then 
do i=l,len 

a ris : ie, 1) ♦--a (is:ie,i)*s(i,i) 

a ^is : ie, i) a (is : i©/ i ) -a { is : ie, i-rl : len) *s (141 : ien, ij 

enddo 

else 

if (Ier.>-32 or len<=20;- then 

ienl-^len/2 

len2'*~lsn-lenl 

else 

lenl^len./3 

len2"^len-lenl 

endif 

call bh-update (a,k,n, s,is,ie, lenl ) 
a(is:ie,l:lenl) a (is : ie, 1 : lenl ) 

-a (is : ie, lenl + i : len) * s ( ler.l+ 1 : len, : : lenl • 
call bh-update (a (1, lenl+l) , k, n, s (lenl+1 , lenl+1} , is, ie, len2} 
endif 
return 
end 
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1 61 



(5 



nothrdfrS^XUc;KiCSljyat>*lfcXU->K»«B 
nb=(n+iWI<-l)/iblk^S£ 



S10 



-S11 



i*<nb-i|rllU^ 





N 





SI 2 



nbase=(nb-1 )*<ibB(S - 



St7 



'S18 



A(nbase+1 ;aribase+ 
1 :nbase+n)C}S|I^J-fcSaLT 

tp(fibase+1 nbase+iblk) lZ^(l> 



A<nbase'« I:nba5e-Hblk. 
1 :rtbase}SJJJ^ (Z) H 

(j^nhase+l 3Q^P>nbase-Hblk 

SifcT<!:?AU.l-nbase)i 
A(ip(j).1:nbase)^ 



'S13 



-S14 



ACnbase+l :n^bas« 

ip(nbase+l:n-1)(=fi 
Aoru^-5„ 



A(nbase*l J*- 
1.1:nbase; 
»«'®ff'^^7hJU©A 



-SI9 



S20 




AttA5L-& 

A(1aij)A<1:r»jpii)>^A 



•S2I 



S22 



C 
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16/ 



1 71 



S30 



S31 
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17/ 



1 8] 



iei=i-Kfl!)«ia$'7a-5^-^-h-e^LfcEI(-?-(D3) 



9 



nbase2=nbase+iWk,len=(r»-nbase2+numthrd-t )/numthrd 
is2=nbase2+<nothrd-l)*len+l . ie2=niinCn.nbase2+nothrd*len) 



S43 




S41 



len=(nba5e-*-nurTithrd- 1 )/numthrd js=(nothrd-1 )*ten+l , 
le=min<nbase.ftOthrd*len) 
3^D*v^© 1 JfeTC @ 5fcS0lstart=nbase+ 1 □•:/'i7(aien=ibfet'^^>„ 




S48 



S49 



nbase2=nbas€+iblk. len=(n-nbase2+numthrd-D/ 
nunrthrdts2=nbase2+(nothrd-l )*len+1. 
ie2-min(n.nbase2+nothrd*ten) 




c 



3 
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[lai 9] 



^3— Ka)iJa3i$3?p-5=-\'-hT*aLfcia(-?-a)4) 



c 



S60 



M 



len=(nbase+numthrd-1 Vnumthrd 
isl=(nothrd~l)*len+l. iel =mln(nbase.nothrd*lenX iof=nbase 



a(isl:ie1,iof+1 :iofH-ib(k)= 
aOsI :ie1 .iof+1 :iof4-iblk)*TRU-U(a(iof4-1 :iof+iblkJof^-1 liof+iblk)) -1 



^ return^ 



S65 



^ it^^U— 5=->d-update ^ 







Ien=(nbase+numthrd-l Vnumthrd 
isl =(nothrd-l )*len+l . ie1=min(nbase.nothrd*len), iof=nbase 












a(iof+1 :iof+iblk,is1 :ie1 )= 
TRL(a(iof+l:iof+iblk,iof+1:iof^-lblk» -1 ♦ 
a(iof»-1:iof4-iblk.is1:ie1) 
TRL Ii]E>5fT5»Ja>T=ftff9a«P» 





^ return ^ 
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2 0] 



S70 



\1 



^xu^^KT'^fif ^2.:ft7cga>jiSjS. *i^.^tt#"r5> 

nbase2=nbase+iblk. Ien=(n-nbase2+numthrd-1 )/numthrd 
is2=nbase2+(nothrd-l )*len+l , ie2=min(n.nbase2+nothrdi'len),iof=nbase 



S71 



a(1:iofas2:ie2)=a(l:iof,iofi-1:iof+iblk) * aOof*-! :io1+iblk.is2:ie2) 



^ return^ 
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[1112 1] 



S75 



len=(nbase+numthrd-1 )/numthrd 
is2=(nothrd-1 )*len+1 . ie2=min(nbase.nothrd*len).iof=nbase 



S76^ 



a(1:iof,is2:ie2)= a(1 :iof,is2:ie2)-a(1 :iof jof+1 :iof+iblk) * a(iof+l :iof+iblk,is2:ie2) 



^return ^ 
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2 2] 



{tin — KG) iias^^? □ — — h-e^ Lfc S ( ^ 0 7 ) 



580-^ 








^ len=(nbase+numthrcl-1)/numthrd 

is2=(nothrd-1 )*len+1 . ie2=(nba&e,nothrd*len).iof=nbase 












a(iof+1:n.is2:ie2)= a(iof+1:n.is2:ie2)-a{ioft-1:n,iof+1:iof+iblk) * a(iof+l:iof+iblkis2:ie2) 







^return ^ 
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V : 22/ 




s(1 :K1 J+l :iblk)=s(1 :iblk)-s(1 J>*s(i J+1 :iblk) 

i=i+1 



^ return^ 
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102 4] 



1t:?';U— ^>e2-update 




tmp=1.0/s(i.i) 
s(i.1:i-1)=tmp*s(i,1:i-l). 

s(r+1 :iblk.1 :t-)=s(i+1 :ibfk.l :i-)-s(i.l :i-1)*s(i+1 :iblk,i) 

s(i+1 :iblkj)='-a(l+1 :iblk.i)*tmp 

a(i,i)=tmp 

i=i+l 



^ return^ 
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2 51 



^13— Ka>42iS$37P— 5^-^— hT?aLfcEl(-€-a)10) 



S95 



396 




S97 





N 


s(1:i-1.1:i-l)=s(1:i-1.1 


:i-1)-s(1:i-l.i)*sai:i-l) 


s(1:i-l,i)- s(1:Hl)*s(u) 


i=l+l 





^ return^ 
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2 6] 



S100 




SI 04 



SI 02 




SI 03 



j5=istartje=istart-1 ♦len-l 
a(js:jle,iS3e)-a(js:Je,isie>-s(is je;i)*a(istart-^i- 1 isJe) 




js=istartje-istart*lefil -1 . js2=Js-nbase je2=js2+lenf -1 
js3'js2+lenl jc3=js3+ten2-l Js4=ista(^len1 ja4=Js4-Hen2-l 
aO&9e.isie}=dUs3ej&:ie>-sOs2ye2Js3-je3)><>a(js4Je4:len.is:ie> 
istart2=isUrt+len1 
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26/ 



2 7] 





S1 15 




is=istartje=istart+lenl - I js2=js-nbase je2=js2+lenl - I 
js3=tstart'*-len1 je3=js2-Hen2- 1 js4=js3-nbaseje4=js4+l8n2- 1 
a(is:iej3js)=a(is:iejs;je)'a(is:iejs3je3>s<j54j5ejs2j82) 
istart2=jstart+lenl 
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1122 8] 



S130- 



barrier f^fflS^'^ 



SI 31 



\J len=(nbase+numthrd-!)/nijmthrcl 

is=(nothrd-l )*len+1. ie=min(nbase,nothrd*len) 
j=nbase+1 




SI 34 




SI 33 



AO.is:ie)^A(ip(j).is:ie)$A*LSX'S> 



i=j+i 



135 



SI 36 



barrier ^Wiit^ 



c 



return 



3 
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m2 9] 



S140- 



barrier ^^^t^ 



SI 41 



len=(n+numthrd-1 )/numthrd 
is=(nothrd-1 }*len+l. ie=min(nbase/iothrd*len) 
j=n-1 



SI 43 



S144 




A(is:iej)i:A(is:ieJpG))$"A*L^X.'g> 



J=j~l 



145 



SI 46 



barrier 1^^$^^ 



( return ^ 
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